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We calculate the production of 2 b-quark pairs in hadron collisions. Sources 
of multiple pairs are multiple interactions and higher order perturbative QCD 
mechanisms. We subsequently investigate the competing effects of multiple b- 
pair production on measurements of CP-violation: i) the increase in event rate 
with multiple b-pair cross sections which may reach values of order 1 barn in 
the presence of multiple interactions and ii) the dilution of b versus b tagging 
efficiency because of the presence of events with 4 5-mesons. The impact of 
multiple .B-meson production is small unless the cross section for producing 
a single pair exceeds 1 mb. We show that even for larger values of the cross 
section the competing effects i) and ii) roughly compensate so that there is no 
loss in the precision with which CP-violating CKM angles can be determined. 
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I. INTRODUCTION 



CP-violation can be accommodated in the Standard Model with 3 families in terms of a 
phase angle in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. There are several initiatives 
to study CP- violation in the S-sector of the CKM matrix. Failure to observe the Standard 
Model predictions implies physics beyond the Standard Model. It is customary to describe 
the .B-sector of the CKM matrix in terms of the 3 angles of the unitary triangle a, f3 and 
7 JO. It is generally expected that a complete determination of the 3 angles will require 
experiments at both S-factories and hadron colliders. 

In this paper we study the measurement of CP-violation angles at the LHC as a function 
of its luminosity. The advantages and disadvantages of the collider have been discussed in 
various publications 0. The LHC is projected to reach a peak luminosity of 10 34 cm -2 s _1 . 
This combined with a large b-quark production cross section at high energy guarantees a 
number of events of the order of 10 10 per year. Even after taking into account the branching 
ratios and detector efficiencies one is still left with a high number of events, of the order of 
10 4 . Therefore the LHC is a S-factory. 

The highest luminosities are, however, reached at a price: the possible presence of mul- 
tiple interactions per beam crossing which could lead to the production of multiple b-quark 
pairs in the detector. The multiple b-pair cross section increases as the number of interac- 
tions per beam crossing and may reach 1 barn at the highest luminosity where it dominates 
the single pair cross section. In extreme cases this effect may increase the number of b- 
pairs produced by two orders of magnitude. Unfortunately, double pair production will also 
introduce an additional source of fake asymmetry. In a double pair event the gold-plated 
tagging e.g. B°B~ — ► J/ipK s l~X as well as the abnormal pairing e.g. B°B + — > J/tpK s l + X 
are possible. Therefore the number of mistagged b's increases since there is no possibility 
of tagging all four 6-hadrons in the event. 

We will in the end conclude that the impact of multiple S-meson production is small 
unless the cross section for producing a single pair exceeds 1 mb. Even in this case, the 
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competing effects of i) a gain in event rates from the additional production of multiple b- 
pair events and ii) the dilution of the CP-signature by the introduction of fake asymmetry 
associated with wrong pairing of b's, are likely to cancel. If one efficiently controls the 
systematics of the mistagging, the sensitivity of the experiments is actually improved. We 
will present results quantifying these statements. 

On the theoretical side, we compute the production of 4 b-quarks in a single pp- 
interaction to leading order in QCD. Its cross section is comparable to the cross section 
for producing 2 pairs of b-quarks by multiple partons when no overlapping events occur. 



a hadron collider event there are 3 sources of events with 2 b-quark pairs: i) higher order 
QCD diagrams, ii) double parton interaction in a single pp collision and, iii) multiple pp 
interactions in the crossing of the two bunches. For high luminosity LHC running mecha- 
nisms ii) and iii) dominate and the cross section for the production of 2 b-quark pairs can 
be written as 



The 3 terms in Eq. [I] correspond, respectively, to the 3 diagrams in Fig.|T]. The first term 
is the cross section for double parton interaction as estimated in || . 02b is the single b-pair 
production cross section. In our calculation we include 02b from perturbative QCD to order 
0(al) ||. R is the radius of the proton occupied by the partons, mostly gluons, producing 
the b-pairs (in our calculation we use R = 1 fm). N is the average number of interactions 
per bunch crossing, 



II. SINGLE VERSUS DOUBLE PAIR PRODUCTION 



At high energies events with more than one heavy quark pair become abundant . In 




(1) 



N = CAta, 



inel- 



(2) 



where C is the luminosity and At the time between collisions. 
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The last term in Eq. [I] is the standard form for the cross section for multiple p-p 
interaction in the bunch crossing with o"j ne ; being the proton-proton inelastic cross section. 
As for the second term in Eq. [I], this is the cross section for the interaction of two partons 
from one proton in one bunch with two protons from the other bunch. It can be understood 
as the interference of the two previous effects. From Eq. |T] we see that 045 increases linearly 
with the luminosity. We present our results on cross sections next. 

In Fig.||].a the single and double bb pair production cross sections in p p collisions are 
shown as a function of \/s. As seen in the figure, at high energy, production by multiple 
parton interaction dominates. 

For the sake comparison we also show the cross section for the the higher-order QCD 
process pp — > bbbb. Its calculation is involved since it receives contributions from 72 Feynman 
diagrams. Our results are in good agreement with the leading-logarithm calculation of 
reference ||. Fig. [| shows several of characteristics of these events. The invariant energy 
of the interaction is around 100 GeV which corresponds to a Bjorken x of about 0.07. Since 
gluon distributions are well tested at this energy we expect our results to be fairly insensitive 
to the choice of parton distributions. The rapidity, transverse energy, and rj-(j) separation 
show that the b pairs are experimentally accessible. However the cross section is two orders 
of magnitude smaller than the single b-pair production cross section. 

As for the multiple parton interaction events, one should expect them to have similar 
distributions as the ones from single events since they are not more than an overlapping of 
those. 

Figure |[b shows the double pair production cross section by multiple parton interactions 
as a function of the LHC luminosity for 2 values of the single pair cross section. The values 
bracket the expectations which strongly depend on the gluon structure function. Multiple 
parton processes may dominate single pair production for £ > 2 — 20 x 10 33 cm~ 2 s -1 for 
o"2fe > 1 — 10 mb. 

It should be pointed out that one cannot expect a to make a reliable estimate of a^b at the 
LHC from leading order perturbation theory. The leading order perturbative calculations, 
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which require inclusion of O(a^) and 0(a^) diagrams, are unreliable because the results are 
sensitive to the assumed quark mass and the renormalization scale. They seem to already 
underestimate the experimentally observed cross section at the Tevatron |7| . Because of its 
intermediate mass the calculation of the bottom quark cross section is believed not to be well 
understood. A perturbative calculation lies beyond the scope of existing QCD technology 
because it requires resummation of large logarithms of 1/x with x ~ m^j ^fs |§. The range 
of values, considered in this paper, is conservative but cannot be guaranteed. 

III. CP VIOLATION MEASUREMENT. RESULTS 

We will illustrate the implications of multiple b-pair production for the study of CP- 
violation using the gold plated channel n\ — > J/ipK s followed by K s — > 7r + 7r~ and J/ip — > 
, which yields information on the angle f3. The b or b nature of the second B is established 
by the sign of the lepton produced in its semileptonic decay. CP-violation results in a non- 
vanishing asymmetry: 

where N +t _ represent the number of events of the type B — > J/ipK s ,B — > J/ipK s , respec- 
tively. The integrated asymmetry is given by 

A=-^BmW), (4) 

where Xd = Am/T ~ 0.69, Am is the mass difference between the two states and T is their 
width. The factor accounts for the possibility of B® oscillation. 

For proton-proton interactions the initial state is not a CP eigenstate, so one expects to 
have a fake contribution to the asymmetry, F, which modifies the asymmetry to 

A ~ a[sin(2/5) + F], (5) 

with a = D-^p. D is related to the so-called "dilution" factor. The equation is valid for 
small (3. F has been estimated to be of order of a few percent of the signal 0. 



We are now ready to analyze the effect of multiple parton interactions on the measure- 
ment of CP-violation at LHC. When including double pair production we must account for 
all possible pairings. We will assume that the trigger is designed to tag a single _B-hadron. In 
the presence of multiple pair production the "other" 5-hadron may, or may not, be assigned 
the "correct" charge. This will increase the fake contribution to the asymmetry. 

Let us define a as the ratio of double to single pair production; 

04& 



J ^ I / AT 1\ . / 1 (N ~ 1) . 



27r.fi: 2 o- iraei V 27r.fi: 2 a ine ; 

where p ~ 0.64 is the probability that two or more B's do not decay semileptonically since 
we are assuming only one prompt lepton in the event. Then the asymmetry can be written 

as: 

_ Nl + aN% -N l _- aN 2 _ 

N\ + aNl + Ni + aNl ' ( ) 

1 2 

where N+_ represent the number of events produced with + or — signature and coming 
from a double (2) or single (1) 6-pair production. We have that 

iV; =1 ' 2 = N(J/ij K s l+X) 1 ' 2 oc {n + +R d h° + R s h s } {1 + A sin 2/3} 



+ { W d h° + W- S h s } {1-^81112/3} 

+ <5 li2 {n+ +R d n° + R s n s } {1 - ^ sin 2/3} 
+ M ^n° + VT/ s n s } {l + T ^sin2/3} 



(8) 



(9) 



iV"i =1,2 = 7V(J/V X s /-X) oc {n" +i? (i n + R s n s } I 1 ~ if% sin 2/3} 

+ { W d n° + W s h s } {1 + ^ sin 2/3} 

+ 5 ii2 {n- + ^n 5 } {1 + ^ sin 2/3} 

+ M WW + WX} {1-^8^2/3} 

The effect of neutral B oscillation in the tagging is parametrized in terms of the W,- L and 
Ri = 1 — Wi coefficients. Wi is the probability of a meson to decay as a B® 

w - - A < io > 



We have used x d = 0.69 and i s = 9 | 

The terms proportional to 5^% account for the possibility of abnormal pairing in the 
double pair processes. We used the same notation as in reference 0: 

2C[N{B + B°)Br{B + -> l+)+ N{A b B°)Br{A b -> /+)] 
2C[N(B~B°)Br(B- -> l~)+ N(A b B°)Br(A b -> Z~)] 
CiV(5 5°)[5r(5 ^/+) + 5r(fi ^/-)] (11) 
2CN{B° s B°)Br(B° s -> /+) 
2CN{B° s B°)Br(B° s -> Z~) 



n 



n 1 
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C[N(B+B°)Br(B + -> Z+)+ N(A b B°)Br(A b -> /+)] 
n - = C[N(B-B°)Br(B- -> Z~)+ N(A b B°)Br(A b -> /")] 

C[N(B°B°)Br(B° -> Z+)+ N(B°B°)Br(B° -f /")] (12) 
CN(B° s B°)Br(B° s -> Z+) 
n s ~ = CN(B° s B°)Br(B° s -> Z~) 
Here n Q and n Q , « = +,—, 0, s, s correspond to "favorable" and "unfavorable" assignments, 
respectively. C is a constant. Even after including all possibilities the asymmetry can still be 
cast in the form of Eq. |5|for small sin(2/3). -F and a must, however, include the contributions 
from multiple pair production. 

Because the initial state is not a CP eigenstate the N(BiBj) are not all equal since 
protons are made of quarks rather than antiquarks and therefore a 6 is more likely to pair 
with a valence quark than a 6. This difference is parametrized |§ in terms of baryon and 
meson fragmentation probabilities l b , l m which are implicitly defined by the equations: 

N(BiBj) oc P{Bi)P{Bj) (13) 



with 



P(B+) = p d (l - Z m ) + 2Z w /3 P(B-) = Pd (l - l b ) 

P(B° d ) = Pd (l - l m ) + Z m /3 P(B° d ) = Pd (l - l b ) 

P(B° S ) = Ps (l- l m ) P(B°) =p s (l- h) 

P{A b )= PA {l-l m ) P{A b ) = PA (l-l b ) + l b 



(14) 



We will assume that pd = 0.38, p s = 0.14 and p\ = 0.1 and take the semileptonic branching 
ratios for A& and B to be equal |J. Clearly for l m = lb = it follows that N(BiBj) = 
N(BjBi) and the fake asymmetry vanishes. The results are not very sensitive to l m but 
depend, in contrast, very strongly on the value of lb- 

In Figs. [| , |5] and ^| we have displayed the value of A(0), the coefficient a, and the fake 
contribution to the asymmetry, F, as a function of the LHC luminosity for different values 
of o"26, (Jinei and the parameters lb and l m . As expected the values of A(0) and F depend 
strongly on the parameter £& while a, which measures the "true" final state asymmetry, is 
rather insensitive to it. As for the effect of multiple pair production, it increases the value 
of the fake contribution to the asymmetry while decreasing the value of the a. As expected, 
it increases the ratio ratio of "fake" to "true" asymmetry. The value of the luminosity at 
which this effect becomes important depends on the value of the For o<ib < 1 mb the 
effect is not very important for any possible value of the luminosity. 

To quantify the impact of double pair-production on the measurement of sin(2/5) we have 
computed the statistical and systematic error on sin(2/3), with 

Vstat - i N 1,10 J 

V ^ * events 



a syst = v /sin 2 (2/3)a2/a 2 + ^, (16) 
where N events is the total number of signal events given by 



(17) 



^events = N(J/lpK s l + ) + N(J /lpK s l~) 

= 2 x N(bb) x P(B d ) x Br(B -> IX) x Br(B d -> J/if)K„) 
xBr(K s -»■ 7T+7T-) x Br(J/^j -> l+l~) x e 
= 6.76 x 10 6 x (£/10 34 ) x (cr 2 &/mb) X (1 + a) events/year 

The e-factor is the product of the triggering, tagging and reconstruction efficiencies com- 
bined with the geometrical acceptance. We use e = 7.7 x 10~ 3 [10|. The systematic errors 



on F and a have been estimated to be of the order of a percent [11]. 
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In summary, the presence of multiple pair production increases the total number of events 
and therefore improves the statistical error. It also increases the fake contribution to the 
asymmetry, F, and, as a result of this, it increases the systematic error. The minimum value 
of sin(2/3) that can be measured with a given statistical significance N a is 



In Fig. [7] we plot sin(2/5) m j„ for 02b = 10 mb as a function of the collider luminosity 
for different values of the systematic errors on F and a for a 3a effect. For comparison 
we show the corresponding value of sin(2/3) min for <j 2 f> = 0.1 mb. As expected the value of 
(5 min and its dependence on luminosity depend critically on the systematic error on the fake 
asymmetry and are less sensitive to other parameters such as the systematic error on a. This 
is easy to understand. Unless the systematic error on F is extremely small, its systematic 
error represents the limiting factor on sin(2/3) m j„ in Eq. [18|. For large values of 026, the 
presence of multiple pair production increases F at higher luminosities and therefore the 
value of the sm(2ft) min . If, on the other hand, one can determine F with good precision, 
the limiting factor in Eq. pi] becomes the statistical error and the increased statistics as a 
result of multiple b-pair production improves the measurement. 



In this paper we have studied CP-violation measurements at hadron colliders. We have 
investigated the effect of multiple b-production as a function of collider luminosity. We 
conclude that multiple b-pair production dominates over single pair production for C > 
2 - 20 x 10 33 cm" 2 s" 1 for a 2b > 1 - 10 mb. 

The presence of multiple pair production produces two competing effects in the measure- 
ment of the CP violating asymmetry. It increases the total number of events and therefore 
it improves the statistics. It also increases the number of mistagged events introducing an 
additional source of fake asymmetry and therefore worsens the systematics. 




(18) 



IV. CONCLUSIONS 
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In the end the impact of multiple _B-meson production is small unless the cross section 
for producing a single pair exceeds 1 mb. In this case the dominant factor weighing the 
competing effects is the error on the determination of the fake contribution to the asymmetry. 
Lower luminosities are advantageous when one cannot measure the fake contribution to 
better than a few percent. However, if a higher precision is achieved, it is still advantageous 
to perform the measurement at a higher luminosity. 
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FIGURES 

FIG. 1. Diagrams for double pair production by multiple parton interaction in high 
energy pp collisions. The three contributions correspond to the three terms in Eq. p] a) 



a 2 Wl a 2 J 1 (JV ~ 1} and rl rr 2 (jV ~ 1} 

a 2b 27TR 2, o) a 2b \j 2nR2 amA , and c) a 2b . 



FIG. 2. a) Single and double bb pair production cross sections in p p collisions as a 
function of ^/s. The shaded band represents the prediction for the single pair production 
cross section a 2b from perturbative QCD to order 0(a^) obtained for several parametrizations 
of proton structure functions and different values of the scale \i and the quark mass m b . The 
upper (lower) edge of the band corresponds to the prediction for the parametrization in Ref. 
|TJ (EHLQ p|]) for the structure functions with m b = 4.5 (5.0) GeV and // = mb/2 (2m b ). 
&AbPERT (dashed lines) is the cross section from higher-order QCD processes pp — > bbbb for 
the EHLQ structure functions with m b = 5 GeV. The upper (lower) curve corresponds to 
ji = Am b (Vs). <JibMP (f uu lines) is the double pair production cross section from multiple 
parton processes (Eq.[TJ). The lower (upper) curve corresponds to N — 1 (50) and the lower 



(upper) values of a 2 b- We have used the prediction of a ine \ from |14[]. Here N is the number 
of interactions per bunch crossing. 

b) Double b-pair production cross section from multiple parton processes (Eq.|l|) at LHC 
{yfs = 16 TeV) as a function of luminosity for (Ti ne i(-\/s = 16 TeV)=85 mb [jnj and At = 25 
ns. The upper (lower) curve corresponds to a 2b = 10 (0.1) mb. 

FIG. 3. Invariant mass (Vs), pseudorapidity (77), transverse energy, and 77 — separation 
(Ar) distributions for the 4-b's events from higher order QCD diagrams. The figures are 
shown for EHLQ structure functions with fi = 4m b . As discuss in the text the distributions 
are not very sensitive to this choice. 
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FIG. 4. Asymmetry for 3111(2/5) = 0, A(0), at LHC as a function of the collider luminosity 
for two values of a 2b . The upper curves correspond to l b = 0. and l m = 0.05. The central 
curves correspond to l b — l m — 0.05. and the lower curves correspond to l b = 0.1 and 
l m = 0.05. In all cases the dark (light) lines correspond to a ine i = 100 (50) mb. 

FIG. 5. Same as previous figure for the fake contribution to the asymmetry, F. 

FIG. 6. Same as previous figure for the coefficient a. 

FIG. 7. Minimum value of sin(2/3) measurable at the LHC with 3a effect as a func- 
tion of the collider luminosity for a 2 6 = 10 mb (solid lines ) and a 2b = 0.1 mb (dotted 
lines), a ine i = 100 mb, and l m — l b — 0.05. From upper to lower the curves correspond to 
a a = a F = 30 , 10 , 1 . 
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